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1.  INTRODUCTION 

The  problem  of  comparing  several  treatments  with  a  control  has  wide 
applications.  Whenever  one  examines  the  possibility  of  replacing  an 
existing  system  by  one  of  several  recently  developed  systems,  the 
central  question  is  whether  the  gain  is  commens'irate  with  the  enormus 
expenditures  that  are  generally  associated  with  such  a  change.  For 
example,  we  may  want  to  consider  replacing  some  equipment,  a  computer  or 
some  drugs.  In  all  such  cases  the  principal  concern  is  whether  the 
existing  system  (control)  results  in  significantly  inferior  performance 
vis-a-vis  some  of  the  new  ones.  A  change  involving  large  financial 
expenditures  may  be  justified  only  when  the  answer  is  truly  in  the 
affirmative.  There  are  plenty  of  such  examples  in  industrial, 
agricultural  and  biological  investigations. 

The  theory  of  optimal  designs,  as  developed  by  Kiefer  and  other 
researchers,  has  traditionally  concentrated  on  developing  good  designs 
for  drawing  inferrences  on  a  set  of  mutually  orthogonal  treatment 
contrasts.  Nothing  much  is  known,  however,  for  the  case  when  the 
contrasts  are  not  mutually  orthogonal  as  in  the  case  of  comparing  the 
control  with  other  treatments.  Till  recently,  there  was  little 
literature  for  this  important  problem  of  comparing  treatments  with  a 
control.  Some  of  the  recent  developments  relevant  to  our  study  are 
Bechhofer  and  Tamhane  (1981),  Constantine  (1983),  Kajumdar  and  Notz 
(1983)  and  Ture  (1982).  Other  notable  references  include  Bechhofer 
(1969),  Bechhofer  and  Nocturne  (1972),  Bechhofer  and  Tamhane  (1983), 
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Bechhofer  and  Turnbull  (1971),  Cox  (1958),  Dunnett  (1955,1964),  Notz  and 
Tatnhane  (1983),  Pesek  (1974)  and  Robson  (1961). 

There  is  a  practical  need  for  a  catalog  of  efficient  designs  for 
coaiparing  two  or  more  test  treatments  with  a  control  under  various 
models  and  settings.  When  there  is  no  nuisance  parameter  in  the  model 
there  is  one  general  easy  solution.  This  is  indicated  in  Section  2. 
Cases  in  which  there  are  one  or  more  nuisance  parameters  in  the  model 

are  very  difficult  to  deal  with.  To  avoid  complications  in  design  and 

its  analysis  later  we  should  run,  if  possible,  the  experiment  under  a 

designed  condition  in  which  the  number  of  nuisance  parameters  have  been 

reduced  to  a  minimum.  Thus  as  a  first  step  in  this  direction  we  need 
efficient  designs  when  there  is  precisely  one  nuisance  parameter  in  the 
model  for  each  response.  This  falls  in  the  category  of  block  designs. 
In  response  to  this  later  need  this  paper  provides  111  efficient 
(formally  difined  as  A-optimal  designs  in  Section  2)  block  designs  for 
cases  in  which  the  block  size  (k) ,  the  number  of  test  treatments  (v)  and 
the  number  of  blocks  (b)  are  in  the  following  ranges 
2  <  k  S  8 

2  i  v  <  10  and  v  -  12,  13,  14,  15,  16,  21,  25 

v  <  b  S  50 

After  formally  introducing  the  concept  of  A-optimal  design  for 
comparing  a  control  with  v  test  treatments  in  Section  2,  we  shall  then 
explain  the  detailed  description  of  our  methods  of  construction.  In 
Section  3,  we  shall  give  a  table  of  our  111  A-optimal  designs. 
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Sufficient  description  is  provided  for  the  construction  of  each  design. 
In  Section  4,  we  shall  summarize  our  findings  and  points  out  some 
interesting  observations. 

2.  A-OPTIMAL  DESIGNS  FOR  COMPARING  TEST  TREATMENTS 
WITH  A  CONTROL 

The  objective  of  an  experiment  is  to  compare  v  treatments  called  test 
treatments  with  a  special  treatment  called  the  control.  Consider  first 
the  situation  where  all  the  experimental  units  are  homogenous  as  far  as 
the  response  to  the  treatments  is  concerned. 

If  y£j  denotes  the  jth  observation  related  to  treatment  i,  then  we 
assume  the  model 

y.j  -  M  +  Ti  +  e£j  (2.1) 

j  "  l,...,n£;  i  “  0,1, ...,v,  0  denoting  the  control.  In  (2.1)  n 

denotes  the  general  mean,  t£  the  effect  of  treatment  i  and  e£j  the 
random  error  which  is  assumed  to  have 

E(e£j)  -  0,  V(e£j)  -  er*.  Cov(e£j,  e£iji)  -  0.  (2.2) 

If  we  are  given  n  experimental  units,  then  the  problem  is  to  determine 
an  optimal  allocation  of  the  v+1  treatments  to  these  units.  Since  we 
are  primarily  interested  in  comparing  treatments  with  a  control,  we 
should  look  at  the  contrasts 

T0  T1 » •  •  • »  T0  ~  Tv* 

If  the  best  linear  unbiased  estimator  of  tq  -  T£  is  denoted  by 

A  A 

Tq  -  r  £ ,  then 

v<r0  “  T£)  -  o’  (1/ng  +  l/n£> 
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To  determine  the  optimal  allocation,  one  sensible  criterion  is  to 

A  A 

minimize  ##v  V(tq-t£).  Thus  we  determine  ng,nj,...,nv 

which  minimizes 

^i“l,...,v  U/HO  +  ^ni^ 
subject  to  the  restriction 

n0  +  ^i-1 , . . .v  ni  “  n* 

In  other  words,  we  minimize 

v(n  -  ^i»i , . . . ,V  ni)_1  +  2i«l,...v  ni_1 
over  positive  integers  nj,...,nv  satisfying 


(2.3) 


V  i  H]  + 


+  nv  <  n. 


(2. A) 


0 

If  the  optimal  solution  to  (2.3)  and  (2.4)  is  denoted  by  n£  ,  then 
the  optimal  nQ  is 

n0#  “  n  ~  2i-lt...V  ni*. 

In  case  n  -  r  ^v(  Vv  +  1),  where  r  is  an  integer,  one  can  show,  using 
some  calculus  that  the  optimal  allocation  is  given  by 

nj#  -  n2#  “  ...  “  nv#  -  r,  ng*  “  r  v'v.  (2.5) 

A 

The  allocation  (2.5)  can  also  be  shown  to  minimize  Max  j<i<v  V (rg  - 
rj) ,  the  maximum  variance  of  the  estimators  of  the  contrasts 

T0-Tl» • • •  *T0"Tv 


Now  suppose  that  all  the  experimental  units  are  not  homogenous  to 
start  with,  but  can  be  grouped  into  b  sets  (called  blocks)  of  k  units 
each.  We  suppose 

k  <  v  +  1  (2.6) 

which  means  that  we  are  in  the  incomplete  block  design  situation  since 
including  the  control,  there  are  v+1  treatments.  Here  the  model  is 
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Yij/  “M+Ti+@j+  eij,  (2.7) 

i  »  0,1,. ..,v;  j  -  1, . . . ,b;  /  ■  l,...,n£j;  n^j  ■  0,1,2,....  There 
is  no  observation  Yjj /  if  n£j“0.  3j  is  the  block  effect  of  the 
jth  block,  and  all  the  other  symbols  have  meanings  similar  to  those  in 
equation  (2.1).  Let  us  denote 

P  “  (  T0"T1’ . r0“rv)' 

the  vector  of  control  -  test  treatments  contrasts.  Our  main  objective 
is  to  allocate  the  treatments  0,1,..., v  to  the  blocks  in  a  way  that 

A 

allows  the  best  possible  inference  on  \p.  In  other  words,  if  p  is  the 
best  linear  unbiased  estimator,  then  we  want  to  'minimize'  the  variance- 

A 

covariance  matrix  V(^),  in  some  sense.  Among  the  three  most  standard 

A 

methods  of  minimizing  V(^),  viz.,  E-optimality  which  minimizes  the 

A. 

maximum  eigenvalue  of  V(^);  D-optimality,  which  minimizes  the 

A  A 

determinant  of  Vty);  and  A-optimality  which  minimizes  the  trace  of  V(\t), 
A-optimality  is  the  only  statistically  meaningful  method  in  this 
situation.  In  this  paper  we  shall  concentrate  on  designs  which  are  A- 

A 

optimal,  i.e.,  which  minimize  trace  V(^). 

Though  some  design  specialists  have  talked  of  block  designs  for 
control  -  test  treatment  comparisons,  a  systematic  study  have  probably 
only  begun  since  the  paper  by  Bechhofer  and  Tamhane  (1981).  The  reader 
is  also  referred  to  their  paper  for  a  survey  of  the  literature. 
Bechhofer  and  Tamhane  (1981)  defined  a  class  of  designs,  known  as 
Balanced  Test  treatment  Incomplete  Block  (BTIB)  designs,  and  discussed 
some  optimal  properties  of  these  designs  for  setting  simultaneous 
confidence  bounds  for  the  elements  of  ip.  A  BTIB  design  is  an  incomplete 
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block  design  in  which  each  test  treatment  appears  in  the  same  block 
which  the  control  the  same  number  of  times  (“Xq)  and  any  pair  of  test 
treatments  appear  together  in  the  same  block  the  same  number  of  times 
(**Xj)  .  Using  the  notation  n^j  introdruced  along  with  (2.7)  we  may 
also  define  a  BTIB  design  by  the  following  relations. 

2j-i . b  nojnij  *  2j«i,...,b  nojnij  “  for 

i  ■  1, . . .  ,v 

2i-i,...,b  nijni'j  “  2j-i,...,b  nijn2j  ■  Msay>»  for 

ii*i  '  “  1 , . . .  ,V 

Note  that  the  usual  BIB  design  is  a  BTIB  design  with  n£j“0,l  and 
X0“X1* 

The  problem  of  determining  optimal  designs  was  considered  by  Majumdar 
and  Notz  (1983).  It  was  shown  that  if  a  BTIB  design  with  the  following 
properties  exist,  then  is  is  A~optimal. 

(i)  n^j  *  0  or  1,  i  *  1,  ...  ,  v,  j  *  1,  ...  ,b 

2j-l,...,b  nij  "  2j-l,...,b  nlj  *  1  “  1 . v 

(ii)  n0£  ”  ...  *  non  ■  L  +  1,  ngN+1  -  ...  -  nQb  -  L 
where  a"L  and  e“M  minimizes 

g(a,e)  •  (v-1) 1 [bvk(k-l)-(ba+e) (kv-v+k)+ (ba2+2ae+e)]-^ 

+  [k(ba+e)~ (ba*+2ae+e)] (2.8) 
among  integers  a  “  0,...t[k/2],  e  "  0,...,b-l,  Ek/2]  being  the  largest 

integer  not  greater  than  k/2.  Note  that  for  each  v,  b  and  k,  there  is 
obviously  an  A-optimal  design,  but  this  result  covers  only  those  cases 
for  which  a  BTIB  design  with  properties  (i)  and  (ii)  exist.  In 
addition,  it  was  also  show  that  '  .s  A-optimal  design  minimizes 


i 

j 
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A  A 

Maxl<i<v  v CT0~Ti) • 

A  BTIB  design  with  properties  (i)  and  (ii)  can  essentially  be  of  only 
two  types: 

Cl)  Rectangular  type. 

1 . b 

1 

L 

L+l 

k 

This  arises  when  M"0.  Here  each  block  has  L  controls.  If  we 
consider  columns  as  blocks,  then  the  designs  looks  as  above,  where  do 
is  a  binary  design  in  b  blocks  of  size  (k~L)  each,  in  test  treatments 
l,...,v.  From  (i)  and  (ii)  it  follows  that  do  is  a  BIB  design.  Here 
'rectangular'  refers  to  do.  A  design  of  this  type  will  be  referred  to 
as  a  R  type,  for  brevity. 

(2)  Step  type. 


CONTROL 


do 


I 


L+l 

L+2 


e 

k 


1 . M  M+l  . b 

1 

• 

• 

L 

L+l 


k 


CONTROL 

_ r 

di 


This  arises  when  M  >  0,  di  is  a  design  in  b  blocks  in  the  test 
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treatments  1,2,..., v.  M  blocks  in  di  have  (k-L-1)  units  each  while  b- 
M  blocks  have  k~L  units  each.  'Step'  refers  obviously  to  di.  One  may 
further  decompose  di  *duUdn,  or  pictorially, 


] 

di  * 

dx  i 

di  2 

dn  and  dn  each  heing  binary  and  proper.  Using  (i)  and  (ii),  it  is 

easily  seen  that  the  test  treatments  in  dn  are  equirepl icated  and  so 

are  in  dn.  Moreover,  if  denotes  the  set  of  blocks  in  d^ 

(i“l,2),  then  Xj  of  (2.7)  can  be  expressed  as 

Xj  “  2  n-nj.-  +  2  n^-n^:  .  (2.9) 

j  «Bi  j  «  B i 

the  same  for  each  i  and  i'  with  i+: '  “  l,...,v.  Also,  Xg  of  (2.7)  is 
a  given  by 

XQ  -  (L+l)  2  nij  +  L  2  n- •  (2.10) 

j  «Bi  j  «Bi 

an  expression  independent  of  i  *  l,...,v.  A  design  of  the  type  will  be 
referred  to  us  as  S  type  for  brevity. 

Let  us  define  r  »  2^^  njj  and  8  *  2j«dn  nlj  *n  a  ^  type 
design.  It  is  easy  to  see  that  for  the  existence  of  a  BTIB  design 
satisfying  (i)  and  (ii)  ,  a  necessary  condition  is  that  there  exists 
nonnegative  integers  r,  <5  and  X£  for  which  the  following  conditions 

hold: 


3 

i 


* 
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vr  -  b(k\  L)  -  M 

vo  -  M(k  -  L  -  1)  (2.11) 

(v  -  1)X2  -  5(k  -  L  -  2)  +  (r  -  5)  (k  -  L  -  1) . 

When  we  are  looking  at  an  R  type  design,  M  ■  0  and  hence  6  is  always  0. 

In  search  of  an  A-optimal  design  with  parameters  v,  b  and  k,  we  first 
obtain  L  and  N  by  minimizing  g(a,e)  in  equation  (2.8).  The  next  step  is 
to  verify  the  necessary  conditions  (2.11).  If  they  are  not  satisfied 
then  there  is  no  BTIB  design  according  to  our  specifications  and  the 
method  fails  to  produce  an  A-optimal  design.  On  the  other  hand,  if  the 
necessary  conditions  are  satisfied  then  there  is  hope  that  a  BTIB  design 
exists,  and  this  will  be  A-optimal  for  the  set  {v,b,k}  we  started  with. 
There  is  no  known  method  which  settles  the  question  of  the  existence  and 

eventual  construction  (if  it  exists)  of  a  BTIB  design  for  a  general  set 

of  paramenters.  So  usually  one  studies  each  case  individually. 
However,  sometimes  the  parameters  {v,b,k,L,K}  allow  a  quick  solution. 
For  example,  if  M~0,  then  we  are  looking  for  an  R  type  design  for  which 
d«  is  a  BIB  design.  In  this  case,  an  appeal  to  some  known  result 
often  settles  the  question.  Usually,  the  BIB  design  may  be  obtained 
from  tables  listed  in  Fisher  and  Yates  (1963),  Rao  (1961),  Hall  (1967) 
and  others.  It  is  interesting  to  note  that  the  d»  part  of  a  S  design 
is  a  design  of  index  \2  38  defined  by  Bose  and  Shrikhande  (1959), 
with  two  different  block  sizes  k~L  and  k-L~l.  In  general,  one  can  think 
of  several  ways  of  constructing  an  S  design,  but  since  they  depend 
heavily  on  the  parameters  v,  b,  k,  L,  M,  we  shall  discuss  the  step 

designs  we  come  across  individually  in  the  next  section. 
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Let  us  now  give  two  examples  of  A-optiaal  designs,  one  each  of  R  and 
S  types: 

Example  1  When  b*14,  v*7,  k“4;  L**l,  M“0.  Denoting  test  treatments 
by  a,  b,  c,  d,  e,  f,  g  and  the  control  by  the  following  R  type  design 
is  A-optimal: 


abcdefgabcdefg 

bcdefgabcdefga 

deffabcdeffabc 

Example  2  When  b*18,  v-6,  k“5;  L“l,  M“6.  Denoting  test  treatments 
by  a,  b,  c,  d,  e,  f  and  the  control  by  -,  the  following  S  type  design  is 
optimal : 


------aaaaaaaabbbb 

aaabccbbbbcdcccccd 

bdbddececddeedddee 

cefeffdfefefffefff 

It  may  be  noticed  in  the  tables  in  the  next  section  that  a  majority 
of  the  A-optimal  designs  is  R  type  with  exactly  one  control  in  each 
block.  Constantine  (1983)  showed  that  in  the  class  of  block  design  in  b 
blocks  of  size  k  each  and  v  test  treatments,  containing  exactly  one 
control  in  each  block,  the  R  type  designs,  which  are  BIB  designs 
augmented  by  a  control  in  each  block  is  A-optimal.  Even  though  many  of 
these  designs  are  A-optimal  in  the  wider  class  of  all  block  designs  with 
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parameters  v,  b,  k,  there  are  instances  where  such  designs  fail  to  be  A- 
optimal.  Consider  for  example,  v"9,  b“48,  and  k*7.  The  design  number 
104  given  in  the  next  section  is  an  A-optimal  design.  However,  a  BIB 
design  with  v»9,  b“48,  and  k“6  augmented  by  a  control  in  each  block  is 

not  A-optimal  since  it  was  not  given  by  our  method  which  produces  all  A- 
optimal  binary  BTIB  designs  whenever  it  produces  one.  Another 
counterexample  can  be  found  in  Constantine  (1983). 

Notz  and  Tamhane  (1983)  and  Ture  (1982)  have  considered  the  problem 
of  constructing  BTIB  designs  which  are  binary  in  test  treatments.  Ture 
(1982)  also  looks  at  designs  which  are  'approximately'  A-optimal. 

Before  concluding  this  section  we  would  like  to  remark  that  the 
Fisher's  inequality,  viz., 
b  >  v+1, 

which  is  valid  for  a  BIB  design  based  on  b  blocks  and  v+1  treatments  is 
not  necessarily  true  for  an  A-optimal  design  in  this  setting.  One  can, 
however,  show  quite  easily,  using  the  standard  technique  of  computing 
the  rank  of  the  incidence  matrix,  that 
b  i  v. 

This  inequality  cannot  be  improved,  in  general,  since  a  symmetric  BIB 
design  augmented  by  a  control  in  each  block  can  be  A-optimal,  as  when 
v“b"k“3 ,  in  design  number  32  in  the  next  section.  However,  for  a  S  type 
design  one  can  show  that  b  -  v+1. 


i'll  ismtanii 
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3.  A  CATALOG  OF  A-OPTIKAL  DESIGNS 
FOR  COMPARING  TEST  TREATMENTS  WITH  A  CONTROL 

Using  the  technique  pointed  out  in  Section  2,  we  searched  for  the 
parameters  of  BTIB  designs  which  made  them  optimal  within  the  following 
practical  constraints, 

2  <  k  <  8 
k  <  v  <  30 
V  5  b  ^  50. 

There  are  precisely  111  sets  of  parameters  satisfying  the  sufficiency 
condition  for  A-optimality  and  the  necessity  conditions  for  BTIB 
designs.  We  were  fortunate  to  be  able  to  construct  all  the 
corresponding  BTIB  designs. 

For  each  triplet  (v,b,k)  the  values  of  L  and  M  (see  the  pictorial 
representation  of  the  design  given  in  Section  2)  were  obtained  by 
minimizing  the  function  g(a,e),  given  in  (2.8),  using  a  computer.  The 
computer  then  verified  the  necessary  conditions  in  (2.11). 

Table  3.1  summarizes  111  cases  for  which  these  necessary  conditions 
were  met,  whenever  M“0  the  corresponding  design  will  be  R  type  and  it 
will  be  S  type  otherwise.  This  is  indicated  in  the  table.  In  the  last 
column  of  this  table  we  have  indicated  how  to  obtain  the  d*  or  d» 
portion  of  the  design,  as  the  case  may  be.  For  those  cases  in  which 
di  could  not  be  described  completely  in  the  limited  space  we  give  the 
entire  layout  in  the  Section  3.3.  For  completeness,  we  also  give  the 
layouts  of  d*  which  happen  to  be  BIB  designs  in  Section  3.2. 
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Throughout  this  section  the  symbol  v2k  denote  all  (  blocks  of  size  k 
in  v  test  treatments.  This  is  clearly  a  BIB  design  -  which  will  be 
referred  to  as  a  complete  design. 

3.1  A  Table  of  A-optimal  Designs. 
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TABLE  3.1 

A  Catalogue  of  A-optimal  Designs 


Design 

number 

V 

b 

k 

L 

M 

Type 

do  or  di 
(The  layout  of  the 

part 

test  treatments  in  this 

design) 

1 

2 

2 

2 

1 

0 

R 

1  copy 

of 

:  221 

2 

2 

3 

2 

0 

2 

S 

dii 

is 

1 

copy  of 

2X1 ,  d i 2  is  1 

copy  of  222 

3 

2 

5 

2 

0 

4 

s 

dii 

is 

2 

copies 

of 

22 1 ,  dii  is 

1  copy  of 

222 

4 

2 

7 

2 

0 

6 

s 

dii 

is 

3 

copies 

of 

221,  dia  is 

1  copy  of 

222 

5 

2 

9 

2 

0 

8 

s 

dii 

is 

4 

copies 

of 

22 1 ,  dii  is 

1  copy  of 

222 

6 

2 

10 

2 

0 

8 

s 

dii 

is 

4 

copies 

of 

22 1 ,  d 1 2  is 

2  copies 

of 

222 

7 

2 

12 

2 

0 

10 

s 

di  i 

is 

5 

copies 

of 

221,  di a  is 

2  copies 

of 

222 

8 

2 

14 

2 

0 

12 

s 

dn 

is 

6 

copies 

of 

221,  di 2  is 

2  copies 

of 

222 

9 

2 

16 

2 

0 

14 

s 

dn 

is 

7 

copies 

of 

221,  dn  is 

2  copies 

of 

222 

10 

2 

17 

2 

0 

14 

s 

di  i 

is 

7 

copies 

of 

221,  di  i  is 

3  copies 

of 

222 

11 

2 

19 

2 

0 

16 

s 

di  i 

is 

8 

copies 

of 

221,  dia  is 

3  copies 

of 

222 

12 

2 

21 

2 

0 

18 

s 

dii 

is 

9 

copies 

of 

221 ,  di 2  is 

3  copies 

of 

222 

13 

2 

24 

2 

0 

20 

s 

di  i 

is 

10 

copies 

of 

221,  dia  is 

4  copies 

of 

222 

14 

2 

26 

2 

0 

22 

s 

di  i 

is 

11 

copies 

of 

221,  dia  is 

4  copies 

of 

222 

15 

2 

28 

2 

0 

24 

s 

dn 

is 

12 

copies 

of 

221,  dia  is 

4  copies 

of 

222 

16 

2 

31 

2 

0 

26 

s 

di  i 

is 

13 

copies 

of 

221,  dia  is 

5  copies 

of 

222 

17 

2 

33 

2 

0 

28 

s 

dii 

is 

14 

copies 

of 

221,  dn  is 

5  copies 

of  252 

18 

2 

35 

2 

0 

30 

s 

di  i 

is 

15 

copies 

of 

221,  di 2  is 

5  copies 

of 

222 

19 

2 

36 

2 

0 

30 

s 

di  i 

is 

15 

copies 

of 

221,  dia  is 

6  copies 

of 

222 

20 

2 

38 

2 

0 

32 

s 

di  i 

is 

16 

copies 

of 

221 ,  di a  is 

6  copies 

of 

222 

21 

2 

40 

2 

0 

34 

s 

di  i 

is 

17 

copies 

of 

221,  di a  is 

6  copies 

of 

222 
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22 

2 

42 

2 

0 

36 

s 

dn  is  18  copies  of  221,  dn  is  6  copies 

of  222 

23 

2 

43 

2 

0 

36 

s 

dn  is  18  copies  of  221,  dn  is  7  copies 

of  222 

24 

2 

45 

2 

0 

38 

s 

dn  is  19  copies  of  221,  dn  is  7  copies 

of  222 

25 

2 

47 

2 

0 

40 

s 

dn  is  20  copies  of  221,  dn  is  7  copies 

of  222 

26 

2 

50 

2 

0 

42 

s 

dn  is  21  copies  of  221,  dn  is  8  copies 

of  222 

27 

3 

12 

2 

0 

9 

s 

dn  is  3  copies  of  321,  dn  is  1  copy  of 

322 

28 

3 

24 

2 

0 

18 

s 

dn  is  6  copies  of  321,  dn  is  2  copies 

of  322 

29 

3 

36 

2 

0 

27 

s 

dn  is  9  copies  of  321,  dn  is  3  copies 

of  322 

30 

3 

48 

2 

0 

36 

s 

dn  is  12  copies  of  321,  dn  is  4  copies 

of  322 

31 

4 

18 

2 

0 

12 

s 

dn  is  3  copies  of  421,  dn  is  1  copy  of 

422 

32 

3 

3 

3 

1 

0 

R 

1  copy  of  322 

33 

3 

6 

3 

1 

0 

R 

2  copies  of  322 

34 

3 

9 

3 

1 

0 

R 

3  copies  of  322 

35 

3 

12 

3 

1 

0 

R 

4  copies  of  322 

36 

3 

15 

3 

1 

0 

R 

5  copies  of  322 

37 

3 

18 

3 

1 

0 

R 

6  copies  of  322 

38 

3 

21 

3 

1 

0 

R 

7  copies  of  322 

39 

3 

24 

3 

1 

0 

R 

8  copies  of  322 

40 

3 

27 

3 

1 

0 

R 

9  copies  of  322 

41 

3 

30 

3 

1 

0 

R 

10  copies  of  322 

42 

3 

33 

3 

1 

0 

R 

11  copies  of  322 

43 

3 

36 

3 

1 

0 

R 

12  copies  of  322 

44 

3 

39 

3 

1 

0 

R 

13  copies  of  322 

45 

3 

42 

3 

1 

0 

R 

14  copies  of  322 

46 

3 

45 

3 

1 

0 

R 

15  copies  of  322 
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47 

3 

48 

3 

1 

0 

R 

16  copies  of  3X2 

48 

4 

6 

3 

1 

0 

R 

1  copy  of  4X2 

49 

4 

12 

3 

1 

0 

R 

2  copies  of  4X2 

50 

4 

18 

3 

1 

0 

R 

3  copies  of  4X2 

51 

4 

24 

3 

1 

0 

R 

4  copies  of  4X2 

52 

4 

30 

3 

1 

0 

R 

5  copies  of  4X2 

53 

4 

36 

3 

1 

0 

R 

6  copies  of  4X2 

54 

4 

42 

3 

1 

0 

R 

7  copies  of  4X2 

55 

4 

48 

3 

1 

0 

R 

8  copies  of  4X2 

56 

5 

10 

3 

1 

0 

R 

1  copy  of  5X2 

57 

6 

37 

3 

0 

33 

S 

SI  in  Sec.  3.3 

58 

8 

40 

3 

0 

32 

S 

S2  in  Sec.  3.3 

59 

9 

24 

3 

0 

18 

S 

S3  in  Sec.  3.3 

60 

4 

4 

4 

1 

0 

R 

1  copy  of  4X3 

61 

4 

8 

4 

1 

0 

R 

2  copies  of  4X3 

62 

4 

12 

4 

1 

0 

R 

3  copies  of  4X3 

63 

4 

16 

4 

1 

0 

R 

4  copies  of  4X3 

64 

4 

20 

4 

1 

0 

R 

5  copies  of  4X3 

65 

5 

10 

4 

1 

0 

R 

1  copy  of  5X3 

66 

5 

20 

4 

1 

0 

R 

2  copies  of  5X3 

67 

5 

30 

4 

1 

0 

R 

3  copies  of  5X3 

68 

5 

40 

4 

1 

0 

R 

4  copies  of  5X3 

69 

5 

50 

4 

1 

0 

R 

5  copies  of  5X3 

70 

6 

10 

4 

1 

0 

R 

R1  in  Sec.  3.2 

71 

6 

20 

4 

1 

0 

R 

1  copy  of  6X3 
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6 

30 

4 

1 

0 

R 

Union  of  the  do 

parts  of 

:  Designs  70  & 

71 

in  < 

this  tabl 

6 

40 

4 

1 

0 

R 

2  copies  of  6E3 

6 

50 

4 

1 

0 

R 

Union  of  the  d« 

parts  of 

:  Designs  70  & 

73 

in  1 

this  tabl 

7 

7 

4 

1 

0 

R 

R2  in  Sec.  3.2 

7 

14 

4 

1 

0 

R 

2  copies  of  the 

do 

part 

of 

Design 

75 

in 

this 

table 

7 

21 

4 

1 

0 

R 

3  copies  of  the 

do 

part 

of 

Design 

75 

in 

this 

table 

7 

28 

4 

1 

0 

R 

4  copies  of  the 

do 

part 

of 

Design 

75 

in 

this 

table 

7 

35 

4 

1 

0 

R 

5  copies  of  the 

do 

part 

of 

Design 

75 

in 

this 

table 

7 

42 

4 

1 

0 

R 

6  copies  of  the 

do 

part 

of 

Design 

75 

in 

this 

table 

7 

49 

4 

1 

0 

R 

7  copies  of  the 

do 

part 

of 

Design 

75 

in 

this 

table 

9 

12 

4 

1 

0 

R 

R3  in  Sec.  3.2 

9 

24 

4 

1 

0 

R 

2  copies  of  the 

do 

part 

of 

Design 

82 

in 

this 

table 

9 

36 

4 

1 

0 

R 

3  copies  of  the 

do 

part 

of 

Design 

82 

in 

this 

table 

9 

48 

4 

1 

0 

R 

4  copies  of  the 

do 

part 

of 

Design 

82 

in 

this 

table 

6 

18 

5 

1 

6 

S 

S4  in  Sec.  3.3 

6 

36 

5 

1 

12 

S 

S5  in  Sec.  3.3 

7 

35 

5 

1 

7 

S 

dn  is  R2  in  Sec.  ! 

3.2,  dn 

is  complement 

of  ' 

4  copies 

9 

18 

5 

1 

0 

R 

R4  in  Sec.  3.2 

9 

36 

5 

1 

0 

R 

2  copies  of  the 

do 

part 

of 

Design 

89 

in 

this 

table 

10 

15 

5 

1 

0 

R 

R5  in  Sec.  3.2 

10 

30 

5 

1 

0 

R 

2  copies  of  the 

do 

part 

of 

Design 

91 

in 

this 

table 

10 

45 

5 

1 

0 

R 

3  copies  of  the 

do 

part 

of 

Design 

91 

in 

this 

table 

12 

33 

5 

1 

0 

R 

R6  in  Sec.  3.2 

13 

13 

5 

1 

0 

R 

R7  in  Sec.  3.2 

13 

26 

5 

1 

0 

R 

2  copies  of  the 

d# 

part 

of 

Design 

95 

in 

this 

table 
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97 

13 

39 

5 

1 

0 

R 

3  copies  of  the 

do 

part 

of 

Design 

95  in  1 

this  table 

98 

16 

20 

5 

1 

0 

R 

R8  in  Sec.  3.2 

99 

16 

40 

5 

1 

0 

R 

2  copies  of  the 

do 

part 

of 

Design 

98  in  1 

this  table 

100 

16 

48 

6 

1 

0 

R 

R9  in  Sec.  3.2 

101 

21 

21 

6 

1 

0 

R 

RIO  in  Sec.  3.2 

102 

21 

42 

6 

1 

0 

R 

2  copies  of  the 

do 

part 

of 

Design 

101  in 

this  table 

103 

25 

30 

6 

1 

0 

R 

Rll  in  Sec.  3.2 

104 

9 

48 

7 

1 

36 

S 

dn  is  2  copies 

of 

R4.  < 

di  2 

is  complement 

of  R3  in  Sec 

105 

14 

35 

7 

1 

14 

S 

S6  in  Sec.  3.3 

106 

15 

16 

7 

1 

6 

S 

S7  in  Sec.  3.3 

107 

8 

28 

8 

2 

0 

R 

complement  of  1 

copy  of 

822 

108 

9 

12 

8 

2 

0 

R 

complement  of  R3  in  Sec 

.  3, 

.2 

109 

9 

24 

8 

2 

0 

R 

2  copies  of  the 

do 

part 

of 

Design 

108  in 

this  table 

110 

9 

36 

8 

2 

0 

R 

3  copies  of  the 

do 

part 

of 

Design 

108  in 

this  table 

111 

9 

48 

8 

2 

0 

R 

4  copies  of  the 

do 

part 

of  Design 

108  in 

this  table 
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3.2  R  Type  Desio£W 

An  R  type  design  is  a  BIB  design  in  the  v  test  treatments  in  b  blocks 
of  size  k~L  each  (called  do)  augmented  by  L  controls  in  each  block. 

In  this  section  we  give  an  exact  layout  of  all  the  BIB  designs  do  used 
in  Table  3.1.  The  test  treatments  are  denoted  by  alphabets.  A  BIB  in  v 
test  treatments  in  b  blocks  of  size  of  k,  each  test  treatment  replicated 
r  times  and  with  each  pair  of  treatments  appearing  in  X  blocks  will  be 
denoted  by  BIB(v,b,r,k,X) .  In  Table  3.1,  These  layouts  are  referred  as 
Rl-Rll.  The  layouts  of  do's  are  as  follows: 

R1  :  BIB(6,  10,  5,  3,  2) 

aaaaabbbcc 

do:bbeddcdede 

efeefdefff 

R2  :  BIB(7,  7,  3,  3,  1) 

a  b  c  d  e  f  g 

do:  b  c  d  e  f  g  a 

d  e  f  g  a  b  c 
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R3 


BIB(9,  12,  4,  3,  1) 


abcadgabcabc 

dtldefbehfdeefd 

ghicfihigigh 

R4  :  BIB(9,  18,  8,  4,  3) 

abcdefghiabcdefghi 

bcdefghiabcdefghia 

d«: 

cdefghiabefghiabcd 

efghiabcdghiabcdef 

R5  :  BIB (10,  15,  6,  4,  2) 

abaaabceeabcabc 

bccbefdfdfdhged 

d»: 

cfgdfhegidgihgf 
•  i  j  h  i  j  h  h  j  j  i  j  i  j  g 
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BIB(12,  33,  11,  4,  3) 

abcdefghijkabedefghijk 
bcdefghijkacdefghijkab 
defghijkabchij  kabcdefg 
hijkabcdefgijkabcdefgh 

abcdefghij  k 
bcdefghijka 
defghijkabc 
11111111111 


BIB(13,  13,  4,  4,  1) 


abcdefghij  kin 
bcdefghijklna 
defghijklnabc 
jklmabcdefghi 


BIB (16,  20,  5,  4,  1) 


d*: 


dhfjhijklmnopegbcabc 
hiikklnnopeegaacddff 
■  nj  plnnopabgdfbhij  k  1 
opoagegabcdcfhijklmn 
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9  :  BIB (16,  48,  15,  5,  4) 

abcdefghijklmnop 
bcdefghi  j  klmnopa 
cdefghijklmnopab 
efghijklmnopabcd 
hijklmnopabcdefg 


abcdefghijklmnop 
bcdefghi  j  klmnopa 
fghijklmnopabcde 
ijklmnopabcdefgh 
klmnopabcdefghij 


abcdefghijklmnop 
bcdefghi  j  klmnopa 
defghijklmnopabc 
hijklmnopabcdefg 
lmnopabcdefghijk 


BIB(21 ,  21,  5,  5,  1) 


defghijklmnopqrstuabc 

ghijklmnopqrituabcdef 

hijklmnopqrstuabcdefg 

mnopqratuabcdefghijkl 

opqratuabcdefghijklmn 
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Rll  :  BIB(25,  30,  6,  5,  1) 


a 

b 

c 

e 

f 

£ 

1 

i 

j 

k 

1 

m 

r 

0 

p 

e 

f 

g 

i 

j 

k 

r 

m 

t 

0 

p 

q 

X 

s 

t 

k 

1 

m 

o 

p 

q 

h 

s 

b 

u 

V 

w 

f 

y 

n 

X 

y 

n 

d 

d 

h 

a 

a 

c 

c 

e 

e 

g 

g 

i 

y 

n 

d 

h 

a 

b 

c 

b 

e 

f 

g 

f 

i 

j 

k 

q 

r 

s 

t 

u 

V 

w 

X 

y 

n 

h 

a 

d 

c 

h 

u 

V 

w 

X 

y 

n 

b 

c 

d 

e 

f 

g 

b 

i 

J* 

i 

j 

d 

1 

h 

a 

0 

p 

q 

r 

s 

t 

u 

V 

w 

j 

k 

k 

m 

m 

0 

P 

q 

r 

s 

t 

u 

V 

w 

X 

1 

m 

1 

o 

P 

q 

r 

s 

t 

u 

V 

w 

X 

y 

n 

3.3  S  Type  Designs 

An  S  type  design  consists  of  a  part  di  which  is  made  up  entirely  of 
the  v  test  treatments  augmented  by  L  or  l+l  controls  in  each  block, 
di  is  the  union  of  dn,  an  arrangement  of  test  treatments  in  blocks 
of  size  k-L-1  and  du  in  blocks  of  size  k-L.  When  k  *  2,  dn  and 
dn  can  only  be  multiples  of  the  complete  design.  But  when  k  >  2, 
construction  of  dn  and  dn  becomes  more  involved.  The  various 
di's  in  this  section  were  constructed  by  different  techniques.  Design 
no.  59  in  Table  3.1  (see  S3  later  in  this  section)  was  obtained  from 
Majuadar  and  Notz  (1983).  In  design  no.  88  in  Table  3.1,  dn  is 

BIB(7,7,3,3, 1)  and  dn  is  BIB(7, 28, 12,3,4) .  In  design  no.  104  in 


Table  3.1,  dn  is  BIB(9,36,20,5, 10)  and  dn  is  BIB (9, 12,8,6,5) .  In 
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design  no.  106  (see  S7  later  in  this  section)  dia  is  the  dual  of  a 
BIB(10, 15,6,4,2) .  In  design  no.  105  (see  S6  in  this  section)  di  was 
constructed  by  removing  the  entire  14  replications  of  one  arbitrarily 
chosen  treatment  from  BIB(15, 35, 14,6,5) .  This  gives  the  required  layout 
in  35  blocks,  14  with  size  5  and  the  rest  with  size  6,  in  14  test 
treatments . 

The  layouts  for  di  for  the  remaining  S  type  designs  (SI,  S2,  S4,  S5 
and  dn  of  S7)  was  constructed  without  relying  on  any  design  available 
in  the  literature.  Let  us  use  the  symbol  A—  (d)  to  denote  the  number 
of  blocks  containing  the  pair  of  test  treatments  i  and  j  in  a  design  d. 
In  di,  A^j(di)  is  a  constant,  independent  of  i  and  j.  Thus,  if  we 
decide  to  construct  du  first  (since  this  is  usually  the  larger  part), 
then  this  can  be  done  quite  quickly  using  several  distinct 
X£j(dia)’s.  However,  this  means  that  there  should  be  many  distinct 
Xij(dn)'s  too;  and  these  are  now  fixed  since  X  - j  (du)  + 

X-(dii)  *  X£j(di),  a  constant.  Constructing  a  du  gets  more 
and  more  difficult  as  the  number  of  distinct  X-j(dn)'s  gets  large. 

Thus  it  is  usually  helpful  to  start  with  a  minimum  number  of  distinct 
X£j(dia)'s  -  usually  2  -  so  that  there  are  also  a  minimum  number  of 
distinct  A—(dn)'s.  The  same  principle  is  helpful  if  we  decide  to 
construct  du  first. 

Denoting  the  test  treatments  by  aplhabets,  we  give  the  layouts  for 
di's  for  k  >  2  as  follows.  They  are  referred  as  S1-S7  in  Table  3.1. 
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SI  : 


a  a  b  c 

dii:  b  d  d  e 

c  e  f  f 


dn  is  2  copies  of  62 2  design  with  the  following  three  additional 
blocks 

a  b  c 
fed 


S2  : 

a  a  a 
di»:  b  d  f 
c  e  g 


b  b  c  c  d 
d  e  e  g  f 
g  h  f  h  h 


dn  is  a  822  design  with  the  following  additional  blocks 
a  b  c  e 
h  f  d  g 


! 
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a  a  b  c  d  e 


di:  aaaabbbbcccddef  fggbf  cdeh 
cdehdeghegi  fifhihi  igfhgi 


aaaaaaaabbbb 


aaabccbbbbcdcccccd 

bdbddececddeedddee 

cefeffdfefefffefff 


aaaaaabbbbcd 


dn:  bbccceccdede 


deddffefffef 


aaaaaaaaaaaaaaaabbbbbbcc 

bbbbbbbbbbccccddcccccddd 

cccccddddeddeeeeddddeeee 


deeffeefffefffffeeffffff 
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a 

h 


f  g  a 
gab 
bed 


1 


g  a  b  c  d  e  f 
£  g  a  b  c  d  e 
d  e  f  g  a  b  c 


c 


g 


n 


j 


d  e  f  g  a 
f  g  a  b  c 
h  i  j  k  1 
i  j  k  1  m 


c  e  f  g  a 
e  £  g  a  b 
n  h  i  j  k 
h  i  j  k  1 


c  d 

d  e 

to  n 

n  h 


k 


j  k 


a  b 


1  n  n  h 


c  a  b  c 


i  j 


a  a 


c 


aabddeefgefgbddd 
bcchifhhhijkcfee 
gfemjgjiillldggf 
kj  inkhkkjnmmhijk 


onmol  lmnooonlnno 
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4.  CONCLUDING  REMARKS 

In  this  paper  we  have  provided  111  A-optimal  designs  for  comparing  a 
control  with  2  to  25  test  treatments  in  blocks  of  sizes  from  2  to  8. 
These  designs  were  obtained  by  using  a  single  method.  Not  surprisingly, 
a  single  method  is  incapable  of  producing  optimal  designs  for  all 
parameter  combinations .  Clearly  there  is  a  great  need  to  develop 
methods  to  cover  the  yawning  gaps  in  the  catalog. 

The  structure  of  A-optimal  designs  present  a  few  unexpected  outcomes, 
rather  contrary  to  conventional  intuition.  For  example,  if  v»k“2  and 
b“10,  one  might  imagine  that  it  is  impossible  to  improve  on  the  design 
consisting  of  5  copies  of 

m 

a  b 

since  our  interest  is  only  in  control-test  treatment  comparisons. 
However  Table  3.1  (design  no.  6)  shows  that  the  A-optimal  design  is 
--------aa 

aaaabbbbbb 

Thus  it  is  an  interesting  problem  to  completely  understand  the  rationale 
for  the  complex  structure  of  these  designs.  One  more  notable  feature  is 
that  the  augmented  BIB  designs  are  often  A-optimal,  but  not  always  so, 
as  pointed  out  in  Section  2. 

For  some  values  of  the  parameter  v.  Table  3.1  contains  no  design  at 


all.  For  these  we  have  to  look  at  larger  values  of  k  and/or  b  to  obtain 
an  optimal  design  by  our  method.  For  example,  when  v*ll,  k-2,  the  first 
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design  was  obtained  when  b*209.  For  the  optical  design  L“0  and  M-99. 
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